IDENTIFICATION DIVISION
PROGRAM-ID ChargeInEMfield
VERSION 20240116
COMMENT A mass m with charge q>0 is located in a static electro-magnetic field [0,0,E], [0,-B,0].
COMMENT Find the trajetory
COMMENT Differential equation (Lorentz force): x'' = omega*z', z'' = a*omega^2 - omega*x'
COMMENT good start settings are omega=0,92, a*omega^2=0,05

ENVIRONMENT DIVISION
ENGINE Anabrid-THAT
TIMEBASE 1ms # corresponds to 1s
REQUIRES INTEGRATOR 4, INVERTER 1, COEFFICIENT 3

DATA DIVISION
OUTPUT OUTPUT.X x
OUTPUT OUTPUT.Z z
COEFFICIENT.1 OMEGA_X # -omega for x'
COEFFICIENT.2 A*OMEGA^2 
COEFFICIENT.3 OMEGA_z # omega for z'

PROGRAM DIVISION
x'' -> INTEGRATOR -> -x'
-x' -> INTEGRATOR -> x
-x' -> COEFFICIENT.OMEGA_X -> -omega*x'
+1  -> COEFFICIENT.A*OMEGA^2 -> a*omega^2
a*omega^2, -omega*x' -> INTEGRATOR -> -z' # input is z''
-z' -> INTEGRATOR -> z
-z' -> INVERTER -> z'
z'  -> COEFFICIENT.OMEGA_Z -> omega*z' = x''

OPERATION DIVISION
MODE REPEAT
OP-TIME 21ms